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1. Introduction 



The kinematic evolution of a car towing n trailers can be described by a Goursat distribution on 
the configuration space M = M. 2 x (S 1 )™+ 1 . A Goursat distribution is a rank-(Z — s) distribution on 
a manifold M of dimension I > 2 + s, such that each element of its flag of Lie squares, 

D = D s C D 3 - 1 = [D S ,D S ] C • • • C D J_1 = [D j ,rP] C • • • C D° = TM 

is of codimension 1 in the following one. 

Since 2000, Goursat distributions were generalized in many works ([7], [12] [13], [14], [15], [20]). 
Special fc-flags (k > 2), which are considered to be extensions of Goursat flags, were defined in 
[7], [14], and [20] in several equivalent ways. All these approaches can be reduced to one transparent 
definition (see [1], [24]). A special fc-flag of length s on a manifold M of dimension (s + l)k + 1 is a 
sequence of distributions 

D s c D s-i = c . . . c l>j-i = [Di,I*>] c • • • C D° = TM 

such that the respective dimensions of D s , D s_1 , ■ ■ • , D° are k + 1, 2k + 1, • • • , (s + l)k + 1, 
for j = 1, • • • , s — 1, the Cauchy-characteristic subdistribution L{D^) of is included in D^ +1 of 
constant corank one, L(D S ) = 0, and there exists a completely integrable subdistribution F C D 1 of 
corank one in D . The integer k is called width. 

The purpose of this work is to show that the problem of modeling car towing n trailers can be 
generalized to the problem of modeling kinematic problem for an "articulated arm" constituted of 
(n + 1) segments in M' £+1 , such that to this model is naturally associated a special fc-flag. 

In the following, an "articulated arm" of length (n + 1) is a series of (n + 1) segments [Mjj Mj+i], 
i = 0, • • ■ ,rt — 1, in E fe+1 , keeping a constant length U, and the articulation occurs at points Mi, for 
i = 1, • • • , n. 

It is proposed to study the kinematic evolution of such a mechanical system under the constraint 
that the velocity of each point Mi is collinear with the segment [Mi, Mj+i], for i = 0, ■ • ■ , n. Such a 
system is also studied in [10] and is called a "n-bar system". In this paper we define precisely the 
kinematic evolution of this mechanical system in term of hyperspherical coordinates and we construct 
the special multi-flag naturally associated to this model. 
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For k — 1, an articulated arm of length n is a modeling problem of a car with n trailers, (see [3]). 
When the number of trailers is large, this problem can be considered as an approximation of the 
"nonholonomic snake" in the plane (see [23] for instance). For k > 1 we can also consider a "snake" 
in R fc+1 (see [22] for a complete description). Again, an articulated arm of length n, for n large, 
can be considered as a discretization of a nonholonomic snake in K fe+1 . For instance, in K 3 , some 
problems of "towed cable" can model in such a way ([17], [23]) 

In section 2, we recall the classic context of the car with n trailers and its interpretation in terms 
of Goursat distribution. The articulated arm system is developed in section 3 and also we show how 
to associate a special multi-flags to such a system in cartesian coordinates. In section 4, we gives a 
version of the kinematic evolution of an articulated arm in terms of angular coordinates and we get 
a generalization of the classical model of the car with n trailers. The last two sections are devoted 
to the proofs of the results. 

2. The car with n trailers 

In this section we will recall some fundamental results about the system of the car with n trailers 
and its relation with the Goursat distribution. All these results are now classical and can be found 
in a large number of papers as [2], [3], [11], [18], [26] and many others . 



2.1. Notations and equations. 

A car with n trailers is a configuration of (n+1) trailers in the R 2 -plane, denoted by M , Mi, • • • , M„, 
and keeping a constant length between each two trailers. It is proposed to study the kinematic 
evolution of the trailer M n with the constraint that the motion is controlled by the evolution of Mo 
which symbolize the car. We will use the same representation as M. Fliess [2] and O.J. Sordalen [26] 
where the car is represented by two driving wheels connected by an axle. It's a kinematic problem 
with non integrable constraints (i.e. a nonholonomic system) due to the rolling without sliding of the 
wheels. The configuration space of the system is characterized by the two dimensional coordinates 
of M n and (n + 1) angles, whereas there are only two inputs, namely one tangential velocity and one 
angular velocity which represent the action on the steering wheel and on the accelerator of the car. 
Consider the system of the car with n trailers and suppose that the distances R r between the different 
trailers are all equal to 1. We choose as a reference point of a body M n _ r the midpoint m r between 
the wheels; its coordinates arc denoted by x r and y r in a given cartesian frame of the plane; 9 r is 
the angle between the main axis of M„_ r and the x-axis of the frame. So, the set of all positions of 
the car with with n trailers is included in a 3(n + l)-dimensional space. This system is submitted to 
2n holonomic links which give, in the previous space, the 2n following equations: 

/ _. \ Xy X 7> — 2. — COS Qy — 2. 

Vr - Vr-1 = sin# r _i 

The configuration space of this problem is a submanifold of dimension (n + 3) which is parame- 
terized by q — (xq, yo,0o, ...,6 n ) where : 

• (xo, ?/o ) ar e the coordinates of the last trailer M n . 

• 9 n is the orientation of the car (the trailer M ) with respect to the x-axis. 

• 9 r , < r < n — 1, is the orientation of the trailer (n — r) with respect to the x-axis. 
The configuration space can thus be identified to R 2 x (S 1 )™" 1 " 1 . 

The velocity parameters are x'o, y'o, #0j • • ' j n . There are only two inputs, namely the "angular 
velocity" w n and the "tangential velocity" v n of the midpoint of the guiding wheels associated to 
the action of the car (see [3] ) . 

Assume that the contacts between the wheels and the ground are pure rolling, it is then submitted 
to the classical nonholonomic links: 



(2) 



x r sin 6 r — y r cos 9 r = 
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There are (n + 1) kinematic equality constraints, one for each trailer. In order to establish these 
constraints, we can represent the points m r , r = 0, • • • , n, in the complex plane, i.e, m r = x r + iy r . 
The geometric constraint between two consecutive trailers is written as: 

m r = m r -i + e l8r - 1 for r ^ 

By induction, we have the following equation: 



r-l 



(3) 



m r = mo + ' 



jo. 



The kinematic constraint of M„_ r is : 



which is equivalent to : 



rh r = X r e 1 ' 



l(e { - M ->m r ) = 

where X(z) denotes the imaginary part of z. Combining this characterization with the derivative of 

(3) and using the linearity of I, we obtain the kinematic constraints: 

(4) — xq sin 9 r + y cos 8 r + T, r ^9j cos(6j — 6 r ) = r = 0, • • • , n 
Combining ni r = \ r e l9r with the derivative of 



m r +i — m r 



1 



we obtain 



and by induction: 



so 



A r = A r+ i(cos6> r+ i — cos# r ) 
A r = A„ cos(6 l „ - 0„_i) • • • cos(6» r+ i - 9 r ) 



rh r = A„( J [ cos(9j - 9j-i))e l6r 

j=r+l 

where A„ = v n is the tangential velocity of the car Mq. 

The evolution of the system of car with n trailers can be given by the following controlled system 
with two controls v n ("tangential velocity") and w n ("normal velocity") of Mq: 

Xq = v cos #o 
y = v sin 9 
9 Q = v\ sin(^i - 6 ) 



(5) 



9 r = v r+ i sm(9 r+ i — 9 r ) 

On-i = v„sin(0„ -9 n -i) 
n = w n 



The "tangential velocity" v r of the body M n _ r is given by : 

n 

v r = | [ cos(9 :j - 9j-l)v n 

j=r+l 
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2.2. Goursat flag. 

Given a smooth distribution Bona manifold M we will use the standard notation [D, D] to denote 
the smooth distribution generated by the vector fields tangent to D and the Lie brackets [X, Y] , of 
any pair (X, Y) of vector fields tangent to D. 

Definition 2.1. A Goursat flag of length s on a manifold M of dimension I > s + 2 is a sequence 
of distributions on M 

D s c D s-i c . . . c d 3 C D 2 C D 1 C D° = TM s>2 (F) 
satisfying the following Goursat conditions 

1) corang D % = i i = 1, 2, • • • , s 

(G) 

2) D 1 - 1 = [D\D l ] i = l,2,---,s 

Each D l {p) is a subspace of T p M of codimension i, for any point p £ M. It follows that D t+1 (p) 
is a hyperplane in D l (p), for any i = 0, 1, • ■ • , s — 1 and p e M. 



Definition 2.2. W^e ca/Z cmy distribution D l of corank i > 2 in a Goursat flag (F) a Goursat 
distribution. 

To each flag (F) of Goursat distributions we associate a flag of "Cauchy-characteristic" subdis- 
tributions 

L{D S ) C LiD 8 - 1 ) C • • • C L{D 3 ) C L(L> 2 ) C L(D 1 ) (L) 

where L(D) is the subdistribution of D generated by the set of vector fields X tangent to D such that 
[X, Y] is tangent to D for all Y tangent to D. L(D) is called the Cauchy-characteristic distribution 
of D. 

Lemma 2.1. (Sandwich lemma)[ll}: Let D be any Goursat distribution of corank s > 2 on a 
manifold M, and p be any point of M . Then 

L(D)(p) c L([D,D])(p) c Dip), 

with dim L(D){p)= dim D{p) - 2, dim L([D, D])(p) = dim D(p) - 1. 

It follows a relation between the Goursat flag and its flag of Cauchy-characteristic subdistributions: 



D s C D 3 - 1 C ••• C D 3 C D 2 C D 1 C D° 

u u u u 

L(D S ) C LID 8 ' 1 ) C L(D S ~ 2 ) C ••• C L(D 2 ) C L{D V ) 

Each inclusion here is a codimension one inclusion of subbundles of the tangent bundle. L(D l ) is 
an involutive regular distribution on M of codimension i + 2. 



2.3. Goursat flag associated to the car with n trailers. 

n 

Let/ r "= J] cosiOj-Oj-!), 

j=r+l 

and v r = /™w„ for r = 1, • ■ • , n — 1 

The motion of the system associated to the car is then characterized by the equation: 

q = w n X*(q) +v n X 2 l (q) 

It is a controlled system with controls v n and w n , (v n is the tangential velocity and w n is the angular 
velocity as we have already seen at the beginning of the section). Each trajectory of the kinematic 
evolution of the car towing n trailers is an integral curve of the 2-distribution, on K 2 x (S 1 )"" 1 " 1 , 
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generated by: 

r = ± 



Xl = cos^o/o n |-+sin^o/o™|-+sin(^i-^o)/r^- + --- + sin(^-^_ 1 )- 



dx UJV d y v " JL de a v 'd6 n . 



1 



The distribution generated by {X^.X^}, naturally associated to the system of the car with n 
trailers, is a Goursat distribution. 



3. Articulated arm 

The purpose of this section is to construct a distribution A, of dimension fc+1, naturally associated 
to an (n+1) articulated arm, and which generates a special fc-flag of length (n+1) on the configuration 
space C = E fe+1 x (S fc )" +1 . Moreover, the kinematic evolution of this arm is an integral curve of A„. 
We begin by recalling the context of special multi-flag in the formalism of [1], [14]. 

3.1. Special multi-flags. 

A special /c-flag of length s is a sequence 

D = D s C D 8 - 1 C • • • C D J ' C ■ ■ ■ C D 1 C D° = TM 

of distributions on a manifold M of dimension (s + \)k + 1 which satisfies the following conditions: 

(i) Di' 1 = [D J ,D j ] 

(ii) D s , D 8 - 1 , • • • , D- 7 ', • • • , D 1 , D" are of respective ranks k + 1, 2k + 1,- • • , sk + 1, (s + l)fc + l. 

(iii) Each Cauchy characteristic subdistribution L(D^) of D- 7 is a subdistribution of constant corank 
one in each Di +1 , for j = 1, ■ • ■ , s — 1, and L(D S ) = 0. 

(iv) there exists a completely integrable subdistribution F a D 1 oi corank one in D 1 . 



Remark 3.1. It should be remarked that the covariant subdistribution F C D is uniquely deter- 
mined by D 1 itself. This covariant subdistribution F is completely described in [1] and [24] where it is 
defined in terms of the annihilating Pfaffian system (D 1 ) 1 - C T*M ([7]). For a complete clarification 
on this fact see [16] . 



Remark 3.2. In the following we mean by a special multi-flag distribution all distribution generating 
a special multi-flag. 

From the definition above, we obtain the following sandwich diagram: 
D s C D 3 - 1 C • • • C D3 C • • • C D 1 C D° = TM 

u u ••• u ••• u 

L{D s - r ) CL(D S - 2 ) C---C L(D J ^ 1 ) C---C F 

All vertical inclusions in this diagram are of codimension one, while all horizontal inclusions are 
of codimension k. The squares built by these inclusions can be perceived as certain sandwiches, i.e 
each "subdiagram" number j indexed by the upper left vertices : 

Di C Di- 1 
U U 
L{D^~ X ) C L(Di- 2 ) 

is called sandwich number j. 
We can read the length s of the special k-fi&g by adding one to the total number of sandwiches in 
the sandwich diagram. 
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Remark 3.3. In a sandwich number j, at each point x G M, in the (k + 1) dimensional vec- 
tor space D^ 1 / L{D^ 1 ){x) we can look for the relative position of the k dimensional subspace 
L(£) J_2 ) I L(D : '~ 1 )(x) and the 1 dimensional subspace / L(D^ 1 )(x): 

either L(D'~ 2 ) / 'L(pi-*)(x) © / L^-^x) = D^ 1 / L(D j ~ 1 )(x) 

or Di/LiDi- 1 )^) C L(D j ~ 2 ) / L(D : >~ 1 )(x) . 
We say that x e M is a regular point if the first situation is true in each sandwich number j, for 
j = 1, ■ • ■ , s. Otherwise x is called a singular point. 

The set of singular points in the context of an articulated arm is studied in [25] and these results 
will be published in a future paper. 

3.2. Special multi-flags and articulated arm. 

The space (R k+1 ) n+2 , will be written as the product Mq +1 x • • • M^ +1 x • • • ffi^+i . Let x, = {x\,--- 

be the canonical coordinates on the space M. k+1 which is equipped with its canonical scalar product 

< , >. (R fc+1 )™+ 2 is equipped with its canonical scalar product too. 

Consider an articulated arm of length (n + 1) denoted by (Mo, • • • M n+ \). In this paper, we assume 
that the distances k are all equal to 1. On (R fe+1 )™+ 2 , consider the vector fields: 

fc+i q 

(6) Z i = Y J (x\ +1 -x\)—toi = Q,---,n 

r— 1 1 

From our previous assumptions (see section 1), the kinematic evolution of the articulated arm is 
described by a controlled system: 

n fe+1 q 

(7) g = yj Mj z 4 +yj Uri+r --_ 

l=0 r=1 OX n+l 

with the following constraints: 

\\xi — x i+ i\ \ = 1 for i = 0- • -n (see [10] chapiter 2). 

Consider the map ^i(x , - ■ ■ , x n+ \) = \\xi — x i+ i\\ 2 — 1. Then, the configuration space C is the 

set 

(8) {(a; , • • • ,x n+1 ), such that ^ l (x , ■ ■ ■ ,x n+1 ) = for i = 0, • • • ,n} 
For i = 0, • • • , n, the vector field: 

is proportional to the gradiant of "Jj. So the tangent space T q C is the subspace of T q (R k+1 ) n+2 which 
is orthogonal to Afi(q) for i = 0, • ■ • , n. 

Denote by £ the distribution generated by the vector fields 

ux n+l ax n+l 

Lemma 3.1. Let A be the distribution on C defined by A(q) = T q CC\£. Then A is a distribution of 
dimension k + 1 generated by 

n n+1 „ 

«+i - II A ^ + 7TT- forr = l---k + l 

i=0 j=i+l UX n+l 

where Aj(q) = - < Afj(q),Afj-i(q) >=< Z j (q),Af j -i(q) > for j = !,-■■ ,n and A n+1 = 1. 
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Proof. Any vector field X tangent to E can be written as: 

n fe+1 



fe+1 Q 



fl/ 

t=0 r=l "+ 1 

On the other hand, on C, a vector fields X is tangent to C, if and only if X is orthogonal to the 
vector fields Wo, • • • ,ftf n - 

For i = 0, • • • , n — 1, each relation < A, A/i >= is reduced to < A^+iZ^+i + \Zi,Mi >= 0, which 
is equivalent to 

(f0) A 4 = X i+ iA i+1 

Similarly, the relation < X,Afi >— induces 

fe+i 



(11) X n =J2^r(x r n+1 -X n , 



r=l 



and from (10) and (11) we get A, = J~J A;A„, for i = 0, • • • n — 1. 

□ 



The properties of A are summarized in the following result, (see also [10] chapiter 2)) 

Theorem 3.1. On C, the distribution A satisfies the following properties: 

(1) A is a distribution of rank k + 1. 

(2) The distribution A is a special k-flag on C of length (n + 1). 



The first part of Theorem 3.1 is a direct consequence of Lemma 3.1. Part (2) will be proved in 
section 6 in terms of hyperspherical coordinates. 



4. The evolution of the articulated arm in a system of angular coordinates 

Given an articulated arm (M , • • • M n+ \) in R fe+1 , we will show that the constraint controlled 
system (7) can be written in the same way as (5) in an adapted system of angular coordinates with 
(fc + 1) controls, namely v n (the "normal" velocity of M n+ \) and the k components of the "tangential 
velocity" of M n+1 (Theorem 4.1). 

4.1. Hyperspherical coordinates. 

The following map 

r(xo, Xl, • , Xi, , X n -|-l) — (,Xq : X\ Xq 7 • , Xi • • • X n +l X n ^j 

implies a global diffeomorphism of (R fc+1 )"+ 2 into itself and T(C) = R fe+1 x (§ k ) n+1 where § fe is 
the canonical sphere in R fc+1 . In this representation, the canonical coordinates on (R fe+1 )™ +2 = 
r((M fc+1 )" +2 ) will be denoted by (xo,z\,--- ,Zi,--- ,z n+ \) so that T is given by x = x and 
Zi = Xi + \ —Xi for i = 0, • • • n. Via this global chart, each point q = (xo, x\, ■ ■ ■ ,Xi, ■ ■ ■ , x n+ i) G C can 
be identified with (xq, Zi, ■ ■ ■ ,Zi,--- , z n -i) £ R fe+1 x (§ fe )" +1 for i = 0, • • ■ , n and C can be identified 
with S = R fe+1 x (§ fe )" +1 . 

We will put on (§ fc )" +1 charts given by hyperspherical coordinates. We first recall some basic facts 
about this type of coordinates. 

The hyperspherical coordinates in R k+1 are given by the relations: 
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(12) 



v = 


p<f> x (9) = p sin 9 1 ■ ■ ■ sin 9 k 


1 sin 9 k 


z 2 = 


p ( j > 2 (9) = psinf? 1 •••sin6» fe - 


1 cos8 k 


< 


p</> 3 (0) = p sin 1 • • • sin 9 k ~ 


2 cos# fe_1 


z k 


pcf> k (9) = psin^cosfl 2 




■ ■ 


= p<t) k+1 {6) = pcosO 1 





with p 2 = (z 1 ) 2 + ■■■ + (z k+1 ) 2 , < 9 k < 2ir and < 9> ' < tt for 1 < j < k - 1. 
We consider $(p,9) = p<&{6) = z, the application from ]0, +oo[x [0, 7r] x ••• x [0, 7r] x [0, 27r] to 
M fc+1 . 

Remark 4.1. T/ie previous expression uses the "geographical" version of hyper spherical coordinates. 
An another version, maybe more usual, can be obtained by taking — — 9 k instead of 9 k and then, 

permuting the functions sine and cosine in each formula. However, our choice is motivated by the 
following fact: the evolution of the articulated arm of length (n + 1) written in a such chart, (see 
(18) ) gives exactly the system (5) for n = 1. 

The jacobian matrix Z?$ of <3E> is: 



fsm.9 1 ■ ■ ■ sin# fe 1 sin# fe pcos^ 1 • • • sin# fc 1 su\9 k ■■■ psinfl 1 • • • sin# fc 1 cos9 k \ 
shift 1 ••• sin cos # fc pcosft 1 ■ ■ •shift' 0-1 cos9 k ■■■ -psinft 1 • • • sin0 fe_1 shift' 



d d d 



\ cosft 1 -psinft 1 / 

fe-i 

It is well known that det(£>$)(p, ft) = (-l)[ fc +V2] ( p )k TJ ( sin gk-iy_ u follows that D $ is invert _ 

i=l 

ible only for < 9 k < 2tt and < ft J < tt for j = 1, • • • , k - 1. 

In the sequence we note {v, O 1 , • • • , fc } the moving frame on 
l>(]0, +oo[x[0, 7r]x- • • x]0, 7r[x]0, 27r[) that is the image, by £>$>, of the canonical frame {-^-, 

Consider a point z = <f>(ft) = <&(1, ft). We note that, in this case, we have 

^ ^ VP dO 1 d6 k > 

where <f> (resp. J^-) is the column vector of components {(f) 1 , • • • , (f> k+1 } 
(resp. {-gfj, • • • , })• 

Note that, in canonical coordinates, this moving frame {u, 9 1 , • • • , fe } can be written 

v ' y J ^ Y dx r y J ^ d9i dx r 

r—1 r—1 

and in hyperpsherical coordinates this moving frame is: 

(14) „(*) = ! and = A 

fc-i 

On the other hand, it is well known that det(£>$)(p, ft) = (-l)[ fc+1 / 2 l (p) fc TJ (sin ft fe ~*)\ It follows 

i=i 

that is invertible only for < 9 k < 2tt and < 9 j < n for j = 1, • • • , fc - 1. 
The column vectors of the jacobian matrix Dl> are pairwise orthogonal and we have 
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IM| 2 = II$H 2 = 1> ||^|| 2 = (sin0 1 ...sin^- 1 ) 2 for. ? =2,...,fc. 
So, at each point z = $(9), the inverse of this jacobian matrix is the transpose of the matrix: 

1 dd> 1 d<t>\ 



(15) 



do 1 



_ .. ll-^-ll dO k I 

\ de i ii — Week W uu / 

This inverse exists only if < 9 k < 2ir and < 9 3 < tt for j = 1, • • • , k — 1. 



Consider a hyperspherical chart y = $(p, 9) + a around some a £ R k+1 and note {v, 6 1 , • • • , fe } 
its associated moving frame. 

Let y' = $'(p',#') + a' be another hyperspherical chart around some other point a' e M fe+1 such 
that its domain intersects the domain of <£. Note {z/, 0' 1 , • • • , 0' fe } the moving frame associated to 
it. So, at each point of the sphere S a of center a of radius 1 in R k+1 contained in the intersection of 
these domains, we can write : 

k 

(16) v' = Av + ^B j e j 

3 = 1 

The components of these vectors are actually the components of the first column vector of the 
matrix o £>$'. According to (15), we get: 

fe+i 

.A(9,6>) = J2<f> r <f>' r 

r=l 

-B\e,9>) = j:%r 



(17) 



•B'{p,Q') = 



1 



fc+i 



2^ ftoj 



\°&\\ ^ 893 
\ee'W r =i 



Remark 4.2. 



the vector ~^^B° Q° is nothing but the orthogonal projection of v' on to TS a . 



4.2. The evolution of the articulated arm on S. 

Coming back to S = R fe+1 x (§ k ) n+1 which is considered as a subset in (R fe+1 )™+ 2 , let § h i = 0, • • • , n, 
be the canonical sphere in K^i 1 . Recall that the canonical coordinates on R^ +1 are denoted by 
Zi = (zj, ■ ■ ■ ,z\,- ■■ , z k+1 ). Given a point a in the sphere Sj, there exists a hyperspherical chart 
Zl+1 = $ 4 (p 4 ,0 4 ) = pi$i(9} , ■ ■ ■ 0f) defined for < 9 k < 2tt and < 9\ < tt, j = 1, • • • , k - 1, 
where $i(0, • • • ,0) = a. So, for a given point q — (xq, z\, ■ ■ ■ ,Zi,--- , z n+ \) 6 S, we get a chart 
(Id — x , $0, • • • i $1, ' • ' ; <&«) centered at q, such that its restriction to pi = 1, i = 0, • ■ • , n, induces 
a chart of 5 (centered at q). For i = 0, • • • , n, in a neighborhood of each Zi + \ 6 R^+i, we consider 
the moving frame 
K i = {v i ,Q},--- ,e k } . 



Remark 4.3. 

(1) Given q — (x , • • • , x n+ \) 6 C, for i = 0, ■ ■ ■ , n, denote by Sj the sphere in R fc+1 of center 
Xi and radius 1. One can put on R k+1 the hyperspherical coordinates yi — <&i(p i} 9i) + Xi. 
As x i+ i belongs to Si, on a neighborhood of x i+ i, we have also the following moving frame 
(again denoted IZi): 

Hi = { Vi ,Ql--- ,o k }. _ 

Note that on x i+ i, the outward normal unit vector ofSi is Vi(x i+ i), and 
{@l(x i+1 ), • • • , Q k (x i+1 )} is a basis T Xi+1 §i. 
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(2) From part (1), if ( = F(q), we get an isomorphism from TcS to 
T X0 R m+1 © T X1 S Q © • • • 8 T Xn+1 § n . 

So we can identify T C S with T xa « n+1 © T X1 S © • • • © T Xn+1 S„. 



Notations 4.1. 

in hyperspherical coordinates, we define on S: 

k+l 

(1) Ai = ^2 4>\-\<fi fori = !,■■■ ,n and A n+1 = 1; 

r=l 

k+l „ 
r— 1 

(3) Z i =J2Bi-^fori = l,---,n 
with: 

' B i ^Y.-Q^tt fori = l,---,n 

r—l * — 1 

k+1 d4> r 

• B l = ,^4-1,1 H ^T 1 ^ fori = l,---,nandj = 2,---,k 

( 4 ) = -qqT> fori = l,---,k,m = 0,---,n 



(5) X^^/;Z !; Km = 0, 



i=0 

m 

with P m = ] [ Aj, for r = 0, • • • , m - 1, and f™ = 1 for m = 0, • • ■ , n. 

j=r+l 

(6) A„ t/ie distribution generated by {X^,X^, ■ ■ ■ ,X%} (with previous notations). 
Remark 4.4. 

(1) For i = 0, • • • , n, consider the column matrix of components (cf>j,--- , and [D^i]^ 1 
the matrix composed by the last k rows of the jacobian matrix of the application 

Finally, denote by [——] (resp. [— ]) the column matrix of components (-^fr,--- i Sk ), 

O&i OX au i 0k> i 

{resp.(-£r, ■■■ , ^t+r)) for i = 0, • • ■ , n. So we can write: 

Zo= f [tH [*o] andZi^ [—— ^[D^]]- 1 ^], fori = 1, • • • , n. 

OX Ovi-l 

(2) According to Remark 4-3 part (2), in the identification : 

T C S = T X0 R n+1 © T X1 S © • • • © T Xn+1 S n , 
the component of X® on T Xi Sj_i is exactly f l n Zi{Q. 

With the previous notations we have the following result 

Theorem 4.1. 

(1) On S, the distribution A„ is the image by T of the distribution A 

where T : C — > S is the diffeomorphism defined at the beginning of Subsection 4-1- 
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(2) The evolution of the articulated arm of length (n + 1) is described in a chart, by the following 
controlled system with k + 1 controls : 



. 1 
X L 

X 2 


= 


j 1 

von 


x k 


fl 


viB\ 




= 


Vl B k 


o\ 




v%+\B\ +1 


k i 




Vi+iB? +1 


o\ 


-1 


= v n B\ 


01 


-1 


= v n B k n 


ol 







n 

where Vi — v n J [ A r and (v$i, ■ ■ ■ ,vgk,v n ) are the (k + 1) controls of the system (18). 

r—i+l 

Moreover, according to Remark 4-3 we have: 

• (vgi,--- ,Vgk) are the "tangential components" of the velocity of M n+ \, namely the 
components, in the canonical basis ofT Xri+ i§ n , of the orthogonal projection of the velocity of 
Af„+i; 

• V{-i is the "normal velocity" of Mi for all i = 1, • • • , n + 1, namely the components of 
the orthogonal projection of the velocity of Mi on the direction generated by Vi-\{xi). 

Remark 4.5. Equations of system (18), fork = 1, are exactly (with the same notations) the classical 
modeling of the car with n trailers ([2], [3], [8], [19], [26]) 

Remark 4.6. According to Remark 3.3 , for k > 2, a point q = (xq, Z\, ■ ■ ■ , z n +i) is singular if and 
only if there exists an index < i < n such that Ai(q) = which is equivalent to [Mj_i,Mj] and 
[Mj,Mj + i] are orthogonal in Mi. In this situation, the velocity all Mj points are zero, for j < i. 
The set of such points is studied in [25] . 

5. Proof of Theorem 4.1 

For part (1), it is sufficient to prove that any germ curve 7 in C is tangent to A at a point geC, 
if and only if, T o 7 is tangent to the distribution generated by X*, • • • , X k } at T(q). 

Let be j(t) — (x (t),xi(t), ■ ■ ■ ,x n+ \(t)) a curve in C defined on ] — e,e[ with 7(0) = q. We have 

(19) Xi = Xi-i + Zi for i = 1 ■ • ■ n + 1 
Assume that 7 is tangent to A. It follows that we have for each t : 

(20) ±i-i = Vi-iZi, for i = 1, • • • , n + 1, 
In view of (20), by differentiation of (19) we get 

(21) ii = Vi-iZi + Zi for % = 1, • • • , n + 1 
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for some Vi E M (which depends of t). 

On the other hand, in R k+1 we also have the orthogonal decomposition 

(22) z i+1 =< z i+1 ,Zi > Zi + Zi i = 0, 
where 

(23) Zi is the othogonal projection of z i+ i on to the hyperplan orthogonal to Zi. 
So we get: 

(24) ±i = v { z l+1 = Vi < z i+1 ,Zi > z l + v t Zi for i = 1, • • • , n 
Comparing (21) and (24), we get : 

(25) Vi-i = Vi < Zi+i,Zi > and z t = v t Zi for i = 1, • • -n 

• for i = n + 1 , we have the orthogonal decomposition 

(26) x n+ i = v n+ iz n+ i + z n+ i 
But we also have 

Xn+l = Z n +\ + X n = V n Z n+ \ + Z n+ l 

So from the unicity of the orthogonal decomposition, we get 

(27) v n = v n+1 and z n+1 = z n+1 

• for i — we have of course 

(28) x = v zi 
From (25) and (27) we get the relation 

n 

(29) r, \ j ■ Zj+i ,Zj > v n for i = • • • , n — 1 

On the other hand, according to Remark 4.3 part (2), as Xi € S,-i, for t = 0, ii can be considered 
as a vector in T x M k+1 and so, in the moving frame IZi, the relations (24), (29) and (26) give rise to 
the following decompositions 

(30) Xi = Vi-\Vi_\{xi) + ViZi for i = 1, • • -n, and i„ + i = u„z/„(a;„ + i) + Z„ +i 
where (rcsp. Z n+ i) belongs to T Xi Sj_i pour i = 1, • ■ • , n (resp. T Xn+1 S n ). 

In fact, according to (23) and Remark 4.3 part (2), Zi is the orthogonal projection 
of (Z)^^) -1 oD$j(i/j(a;j + i)) on T Xi Sj_i, for i = 1, • • • , n . Taking in account Remark 4.2 and (17), 
we get the expression of Zi given in Notations 4.1 for i = 1, • • -n . Moreover, the i e component 
£i of DT(j(0)) is Vi-\Zi in hyperspherical coordinates for i = 1, • • - n . On the other hand, in the 
hypcrspherical $„, the vector Z n+ \ can be written 

Z n+ l — ID, 

Finally, from (29) and the value of Aj, we get: 

Vi = f % n v n for i = 0, • • • n 
On the other hand, from (28), as vector in T Xo R m+1 , we have the following decomposition: 

m 3 

r=l 

So, according to Remark 4.3 part (2) and Remark 4.4 part (2), in hyperspherical coordinates, we 
finally obtain that, DT(j(0)) is tangent A„ at £ = T(q). 



ARTICULATED ARM AND SPECIAL MULTI-FLAGS (CORRECTED VERSION) 



13 



Conversely, let be S(t) = (xo(t), z\(t), • • • , z n+ i(t)) a curve in S, denned on ] — e, e[, and such that 
5(0) = C = r(g). Assume that, in hyperpsherical coordinates, 6 is tangent to A„. and so S satisfies 
the differential equation (18). According to the definitions of Zi in hyperpsherical coordinates and 
taking in account (23) and Remark 4.3 part (2), there exists a curve v n (t) in R such that, for any 
i = 1, • • ■ n, we have 

n 

Zi = [( Yi < z 3+l> z 3 >) v n][Zi+l- < Z i+1 ,Zi > Zi] for i = 1, • • • n 

and 

n 

(31) x = [(JJ < >)u/,]^i. 
It follows that for i = 1, • • • n we have: 

« n 

(32) ij = i + ^ Zj = ( Y[ < z 3+-t> z J > v n)Z l+ l 

j=l l=i+l 

so, according to (31) and (32), we obtain that -f(t) — o S(t) is a curve in C which is tangent 
to A. 

Taking in account part (1) of Therorem 4.1, the kinematic evolution of the articulated arm is a 
controlled system on S which is exactly (18). However, for the completeness of the proof of this 
result, we must prove the interpretation of the control in terms of the components of the velocity of 
M h i = !,■■■ ,n + l. 

From (30), according to Remark 4.3 part (2), in hyperpsherical coordinates $„ around x n+ \ G 
the velocity x n+ \ of M n+1 can be written: 

X n +1 = V n U n (x n+ l) + V" V gj -^-r 
~[ n dV J n 

Where v n is the "normal velocity" of M n+ \ and {v$i, ■ ■ • ,Vgk) are the "tangential components" of 
the velocity of M n+1 . Moreover, according to the value of X®, the contol parameters are exactly 
{v n ,v e i,- ■ ■ ,Vgh). 

In the same way, for any i = I,-- - ,n, in hyperspherical coordinates around Xi G R k+1 the 

velocity ±i of Mi has the decomposition 

±i = Vi-ii>i-i(xi) + v t Zi 

where Zj belongs to T^Sj-i. So, Vi-\ the "normal velocity " of Mi 

6. Proof of Theorem 3.1 

We will see that the distribution A„ generates actually a special fc-flag of length (n + 1) on a 
k(n + 2) + 1 dimensional manifold. Let's introduce the following notations: 

• A m is the distribution generated by {X^, X^, ■ ■ ■ X^} for m = 1, • • • , n; 

• D m+1 is the distribution generated by X^ and {Xj, • • • , X!f m < j < n} 
for m = 0, • • • , n; 

• D° =TM 

• E m+1 is the distribution generated by {Xj, • • • , X k , m < j < n} 
for m = 0, • • • , n. 
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Proposition 6.1. A„ is a special k-flag distribution. More precisely, it satisfies the following 
properties: 

(1) For m = 1, • • • , n + 1, the distributions D m and E m , are of respective constant dimensions 
(n — m + 2)fc + 1 and (n — m + 2)k; 

(2) for to = 1, • • • ,77 + 1, E m is an involutive subdistribution of D m of codimension 1. More- 
over [E m+ \D m+1 ] C D m for to = l,--- ,77. Actually E m+1 is the " Cauchy- characteristic 
distribution" ofD m for m = !,-■■ ,n ([13], [14]); 

(3) [D m+1 ,D m+1 ] = D m for all m = 0, • • • , 77; 
(4) 

A„ = D n+1 c • • • C D m c • • • C D 1 c D° = TM 

u ••• u ••• u 

E n+1 C • • • C E m C • • • C E 1 

Proof of Proposition 6.1. 

It is sufficient to show the property (4). The inclusions [E m+1 , D m+1 ] C D m for 777 = 77, • • • ,0, 
are an easy consequence of (4) and the properties (1), (2) and (3) are always true, according to the 
definition of spaces E m , D m and A m . 

d d 

Denote by An the distribution generated by \ -t-t, • • • , — rrr r- 

1 dxl dx* +1 

For all 777 = 1, • • • , 77 + 1, we have : 

D m = E m+1 e i = D m+1 + 

£)"»+!] contains the space generated by D m+1 and the Lie brackets [X l m , X .], for i = 1, • • ■ , fc. 
We will show by induction that, in fact, they are generating [Z) m+1 , D m+1 ] modulo D m+1 . 

For all m = n, ■ ■ ■ , 0, we have X . = A m X^ n _ l + Z m . It results from the definition of Ai, Zi and 
X° that [X m , X^_A = 0. So we have 

[ X mi X m\ = X m{Am)X m _ l + {X* n , Z m }. 

For j = 1, • ■ • k + 1, consider the vector fields: 



V^' V° I \ " 1 d( Pm-l y 



m—1 

r =i ii ae^"_i II ""m- 1 
If we set $ m _i(/9, 6> m _i) = p$ m -i(9m-i), then we have the relation: 

[^m-l] = {D^ m -l\ 1 [X m -l] 

where the vectors column [F m _i] and [X m _i] have {Y^ l _ 1 , ■ ■ ■ ,Y^\} and 

{X^ l _ 1 ,-- - , X TO-1 } as components respectively. It results that - i^m^i} i s a basis of 

A m _i. 

For 777 = 0, • • • , 77, we note [£)$ m ] the jacobian matrix of $ m (p, 9 m ) = p<& m {9 m )- 

The following decompositions occur: 
fe+i 

3=1 

l x L x l] = E i for all 777 = 1, ... ,77 

By similar way for -D<3? m , we can show that the family of vector fields 

{ X mi \ X rrn X m\i i \ X rm X m\\ 

is also a basis of A m _i . This result is also true for m = 0. 

Since L> m+1 = £ m + 2 © A m , the space [D m+1 , D m+1 ] contains E m+2 , all vectors 

X^, X x ml ■ ■ • , X^ and the Lie brackets [X^, X£j, • • • , [X^,X£J. Also, all the Lie brackets [X^, X£j 
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are zero for r = m + 1, • • • , n and j — 1, ■ • ■ , k since X^ does not depend on variables O^, for 
r = m + 1, • • • , n and j = 1, • ■ • , k. The other Lie brackets [X^X^] are zero for r = m + 1, • • • , n 
and i, j = 1, • • • , k. 

Since {X® nl [X^.X^, • • • , LY^,X^]} is a basis of A m _i, then we have 

[D m+1 ,D m+1 ] = E rn+1 A m _x = D m 
which completes the proof of proposition 6.1 and Theorem 3.1. □ 



Comment 6.1. 

Given two integers p and m such that 1 < p < m < n, we can look for the motion of a "sub-induced 
arm", which consists of segments of the original arm between M p _i and M m+ \ included. We can then 
study the motion of M p _\ as the motion of the extremity of this sub arm for the motion commanded 
by the segment [M m ; M m+ i] . We put h = m — p + 1, and we write H p , m for the canonical projection 
from R k+1 x (S k ) n+1 onto R fe+1 x (§ k ) h+1 defined as U p ^ m (x, Z\, ■■ ■ , z n+1 ) = (x p -\, z p -\, z p , ■ ■ ■ , z m ) 
where x p ^\ are the Cartesian coordinates of M p _\. 



The evolution of the extremity M p _i of this articulated sub-arm , controlled by the movement of 
[M m ; M. m+ i], is a solution of the following differential system (with notations of Theorem 4-1): 



(33) 



' x 1 


= v P -i<t>l-i 


T 2 

x p-l 


= v P -i<t>l-i 


+ fe+l 


= v P -^ k +_\ 


9 p -i 


= v p B p 


ok 

tfp-l 


— v p B k 


crB\ 


= v i+1 Bj +1 


6 k = 


v t +iB k +1 


#m-l 


= v m B x m 


ok 


= v m Bi 


Qln = 


V0l m = v m +i 


ok _ 
^ra - 


v e ^ = v m+ i 



r>k 

n m+l 

It is a controlled system on R fc+1 x (§ k ) h+1 (h = m - p + 1) : 

k 

'q = u Xl + u ^ l h 



with controls uq — v m , and Ui = v S i , for i = 1, 



i=i 



q = Ii Pym (q) 



X°h = E M + ffn'Zp-l et 4_! = ]T 



;=i ox p-i 



We denote by Ah the distribution generated by X® and X^, ■ ■ ■ ,X k n . 



This comment will be used in a future paper about singular sets of special flags and their inter- 
pretations in terms of singularities kinematic evolution of an articulated arm. 
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